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ABSTRACT 


This paper deals with the synthesis of robustly stabilizing compensators for interval 
plants, i.e., plants whose parameters vary within prescribed ranges. Well-known H°° 
methods are used to establish robust stabilizability conditions for a family of plants 
and also to synthesize controllers that would stabilize the whole family. Though 
conservative, these methods give a very simple way to come up with a family of 
robust stabilizers for an interval plant. 

1. INTRODUCTION 

In parametric robust control theory, a plant with uncertain parameters is modeled by 
an interval transfer function, i.e., a transfer function whose numerator and denomina- 
tor polynomial coefficients vary within prescribed ranges. Thus, we have a family of 
plants and the whole family needs to be stabilized with the help of a fixed controller. 
Stability analysis, in the area of parametric robust control, is well - developed. The 
problem of determining whether an interval control system is robustly stable, is solved 
for single input multiple output and multiinput single output systems by an extremal 
theorem that appeared in [1]. Given an interval plant and a controller, this theorem 
reduces the problem of checking the stability of the system to the problem of checking 
the stability of a prescribed set of extremal line segments, the number of which is in 
dependent of the dimension of the uncertainty box in parameter space. 

However, a synthesis synthesis theory for robust stability of interval plants is yet 
to be developed. Here an attempt is made to convert the synthesis problem for an 
interval plant to an H°° synthesis problem which is solvable. The extremal theorem 
is used to convert the parametric uncertainty to an uncertainty band which is used 
in the H°° synthesis method. The method developed in [2] is used. This method 
uses Nevanlinna - Pick theory to derive a condition for robust stabilizability and to 
synthesize a family of robust stabilizers. 

The method developed in this paper, though simple and elegant, is conservative. 
This conservativeness gets introduced while choosing the uncertainty band function 
needed for the H°° methods. However, this conservativeness could possibly be over- 
come in the future by using some loop-shaping concepts or some results of approxi- 
mation theory. 

2. PROBLEM FORMULATION 

We consider a single input single output (SISO) interval control system as shown 
in Figure L The plant P(s) is an interval plant built around the nominal plant 
P 0 (s). Let the entire family of plants be denoted by V. Let the perturbations of the 
coefficients of the plant be parametrized by e, i.e., the perturbations are characterized 
by a box in the parameter space. Each side of this box has a length e, weighted by a 
weighting factor. The values of these weights depend on the scaling inherent in the 
physical parameters. 

The problem of obtaining robust stabilizability conditions and also that of syn- 
thesis of robust stabilizers for interval plants is as yet unsolved. In this paper, the 
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H°° robust stabilizability conditions, as dealt with in the following sections, are used 
to obtain similar conditions for interval plants. Also, the H°° synthesis methods are 
used to design robust stabilizers for interval plants. 

In H°° control theory, a plant with perturbations is characterized by a nominal 
model and an uncertain band function, i.e., the perturbed plant P(s) is said to belong 
to the class C(P 0 (s), r(s)) where it is assumed that 

1) P(s) has the same number of unstable poles as Po(a), 

2) | P(jw) - P 0 (ju>)\ < |r(;u/)|, |r(jw)| >0, Vw. 

So, we find that the H°° uncertainty is modeled by the stable, uncertainty band 
function, r(s), while the parametric uncertainty appears in the coefficients of the 
transfer function model of the plant. If, in any way, we could link up the two, then 
the H°° synthesis methods could also be extrapolated to the parametric robust control 
case. So the essential problem here is to obtain a stable r(s) corresponding to an e. 
Then, we could go ahead and find the maximum e, fmax, for which the plant is robustly 
stabilizable and also find a family of controllers that stabilizes the entire family of 
plants, V, in the e m « - box. 

3. ROBUST STABILIZATION USING NEVANLINNA-PICK THEORY 

In H°° synthesis, the problem of robust stabilization of a class of plants C(Po(s), r(s)) 
is reduced to the problem of finding an SBR function it(s) which interpolates to 

u(ai) = fa, i = 1, . . . Re[a,] > 0, |/?,-| < 1 (1) 

This synthesis method is discussed in [4]. The general closed-loop system as shown 
in Figure 1 is considered. A robust stabilizer C(a) needs to be synthesized for the 


+ 


a 


C(s) 


P(s) 


Figure 1: General Closed-loop System Configuration 


interval plant .P(s). Here, 


u(ai) = 


^m(a«) 

PoK) 


= 0i 


( 2 ) 


where the are the poles of Po(^) in Re[s] > 0. It is assumed here that Po( J ) b as no 
poles on the jco axis and r m (s) is nothing but a minimum-phase H°° function such 
that 

km(jw)| = |r(iu/)| 
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Table 1: Fenyves Array 


ai a 2 

a 3 ... 

Oil 


Pi Pi, 2 

Pm • • • 

Pl,l 

«l(a) 

P 2 

P2,3 

P2,l 

u 2 (s) 



Pi 

«*(*) 


and 

P 0 (s) = B(s)P 0 (s) 

where 

B(s) = TT ( — — - J (Blaschke product). (3) 

AA \a< + s/ 

Now, the problem of finding the SBR function u(s) which interpolates to u(a{) = 
0i is solved by using the Nevanlinna - Pick interpolation technique. This technique 
is discussed in detail in [2]. In this technique, the algorithm for finding such an u(s) 
is as follows: 


1. First the Fenyves array is formed as shown in Table 1. where Ui(s) = u(s), p\ = 
and pij = 0j] for the other rows 


Pi-1, i ~~ Pi-1 

Pi = i _ 

+ l<i<l 

(4) 

1 - Pi-iPi-i,i 

ai - «i-i 


Pi-1, j Pi— 1 

Pi,j — 

+ 1 < i < j < l 

(5) 

1 - Pi-iPi-ij 

a; - «i-i 



2. A Schur function ui(s) is parametrized in terms of an arbitrary Schur function 


ui(s) = 

such that ui(s) interpolates to 


Pi + «i+i(s) 

( 

V*+Ot| 

) 

1 + Pi 1 


Ul +1 (sY 


( 6 ) 


M a i) = Pi 


3. The functions u;_i(s), u/_ 2 (s), .. ., u 2 (s), u-l(s) are computed by repeating the 
linear fractional transformations 


Ui_i(s) 


l+*-i (S§S)«<(*) 


(7) 


A family of u(s), parametrized in terms of an arbitrary Schur function iq + i(s), 
can be obtained this way. 
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A condition to test whether a solution exists is given by the following theorems. 

Theorem 1. The Nevanlinna-Pick problem admits a solution iff the modulus of all 
the elements of the Fenyves array are less than one , i.e 


Ip* I ^ \Pij\ ^ 1- 

From the Maximum Modulus Theorem, if \pk\ = 1 for some k then only one solution 
Ufc(s) = pk exists if pk,j = Pk • /n this case we have 

= \u k - 7 (jv)\ = ■ - • = \u\(ju>)\ = 1 Vw 

Theorem 2. The Nevanlinna-Pick problem admits a solution iff the Pick matrix P, 
whose elements pij are given by 

i -Mi 



is non-negative definite . 


This technique can be modified for plants with poles at the origin and is dealt 
with in [4]. Also, in the case where some points at; are complex, the Nevanlinna - 
Pick algorithm can be modified accordingly and is reported in [5]. 

Now, once 11 ( 5 ) is obtained using Nevanlinna-Pick interpolation we can compute 
q(s) from 

W = ( 9 ) 


and the controller C(s ) from 




C( s ) - ***) 

W l-Po(s)q(sY 


( 10 ) 


To get a proper controller, u(«s) needs to be strictly proper (i.e., u(oo) = 0) when 
the relative degree of r m (s) is 1. If r m («s) has a relative degree greater than one the 
interpolation problem at 00 is more complex, but solvable. 

When Po(*s) is stable, there are no interpolation conditions and it is only neces- 
sary to choose ^( 5 ) as any SBR function. Then q(s ) can be obtained from 


9(3) = 


r m(a) 


( 11 ) 


and the compensator is 


C(s) 


sif) 

1 - Po{s)q(s)' 


( 12 ) 
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4. MAIN RESULT 


First, let’s see how we could find a stable uncertainty function r(s) in the H°° con- 
trol domain that would correspond to an e in the parametric domain. For this, the 
parametric perturbation of the plant is represented as an additive unstructured per- 
turbation, AP(i), i.e., 

A P(s) = P(s) - P 0 (s) (13) 

where P(s) is the interval plant. Obviously, AP(s) will be a family of transfer func- 
tions. Now, we can say from [1] that at each frequency, the maximum magnitude, 
in the family of transfer functions P(s), will correspond to a point on one of the 
extremal segments, denoted as the E - segments of P(s). The same result can be 
applied here to show that the maximum magnitude of A P(a), at each frequency, will 
also correspond to a point on one of the E-segments of P(s). This way, we can search 
over the E - segments of P(s) at all frequencies and get maxp 6 p |AP(, 7 ’uj)| for all u>. 

Now, one way of getting an r(s) is to choose 

(14) 

the maximum uncertainty ball over all frequencies. The problem with the choice of 
constant r(s) is that it gives a conservative bounding of A P(s), thus producing a 
conservative result, which means a lower value of f m »i- The conservativeness can be 
minimized by introducing poles and zeros in r(«) and shaping it such that |r(ju;)| 
approximates ma xp 6 p |AP(ju;)| at each w as closely as possible. 

So, for each e, the maximum magnitude plot of AP(s) can be obtained by 
searching through the E - segments of P(s), the maximum magnitude of A P(s) 
increasing with e at all frequencies, and an r(s) can be loop - shaped such that the 
|r(ju/)| plot approximates the maxp e p |AP(jw)| plot from above. 

Let us choose r(s) to be a constant as shown in eq. (14). Then r can be obtained 
for each e. The plot of r vs e will have the controller design information. 

4.1. Using Nevanlinna - Pick Method 

The Nevanlinna-Pick technique uses this r to check for the existence of a robustly 
stabilizing controller. Then, with the help of Theorems 1 and 2, it is possible to find 
a maximum r, r max such that the H°° problem has a solution for all r < r m ^ T . This 
means that r max is the largest value of r for which there exists a compensator C*(s ) 
which stabilizes all plants in the class C’(Po(s),r). The elements of the Pick matrix 
and Fenyves array being functions of r, this happens iff the Pick matrix is non - 
negative definite or the Fenyves array has all its elements of modulo less than 1 for 
all r < r m „. Therefore, r max is that r for which |p fc | = 1 for some k in the Fenyves 
array. 

So the following algorithm can be used to find 

1) All the elements of the Fenyves array are computed, starting with the first row. 

Obviously, the elements are going to be functions of r. An initial arbitrary choice 

is made for r. 


r(s) = 


max | AP e (ju;)| 
Per v 
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2) If \pk\ — 1, for some row k, go to step 3. If |/>i| < 1 for all i, increase t and go to 

step 1. If \pk | > 1 for some k, decrease r and go to back to step 1. 

3) Let Uk(a) = pk, and compute the other functions «fc_i(s), Ufc_ 2 (s), . . ., «i(s). 

4) A solution to the interpolation problem is then given by 

u(s) = Ui(s) 


with an H°° norm equal to 1. 

Now let P(s) be an interval plant with an e - box perturbation around the nominal 
plant P 0 (s). For each e, the E - segments of P(s) can be searched at each frequency 
to get maxpg'p | AP(jw)| for all u. r is nothing but the maximum of these values over 
all u>, and obviously increases with increasing e. So a maximum e, can be reached 
for which t = r max , i.e., 




max \AP em „{ju)\ 


Hence, by the Nevanlinna - Pick conditions, e max is the largest box in the coefficient 
space for which the interval plant P(s) is stabilizable. Of course, this value of e max is 
a conservative estimate as r(s) is chosen to be a constant. The compensator C*(s ) 
will stabilize the family of plants V € where 


Pc — {Pd + AP e | € < £ni»x}' 


(15) 


The following examples will illustrate the above theory. Let us start with an 
example for a stable nominal plant. 


Example 1. Let the nominal plant be 

p, \ _ a + l 

01 } s 3 + 8s 2 + 22s + 20 

with poles at —2, —3 ± j, all in the LHP. The perturbed plant is 

p( 3 ) = 3 + . a 

K J s 3 + bs 2 + cs + d 

and the perturbation of the coefficients about the nominal is given by 

d 6 [1 — e, 1 + c], 

b G [8 — e, 8 + e], 
c G [22 — e, 22 + e], 
d G [20 — e, 20 + e]. 

Since, it is assumed that the number of unstable poles of the plant should remain 
unchanged, in this case it is required that e < such that the entire family of 
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plants is stable. This initial bound, e m „ on e can be found by checking the Hurwitz 
stability of the Kharitonov polynomials 

K d = s 3 + (8 + e)a 2 + (22 - e)a + (20 - e) 

K J2 = a 3 + (8 + e)a 2 + (22 + e)a + (20 - e) 

Kjs = s 3 + (8 - e)s 2 + (22 - e)s + (20 + e) 

K d 4 = s 3 + (8 - e)a 2 + (22 + e)a + (20 + e) 

From the Routh - Hurwitz table for K dl , K d2 , K t o, K di , it can be shown that e™* = 
6.321, i.e., the entire family of plants will be stable as long as e < 6.321. 

Now for each e, all the E - segments will have to be searched for all ui to get 
maxp e 7 > |AP(jw)|, the H°° norm of which will be chosen as r corresponding to that e. 
Once an r is obtained, a robustly stabilizing compensator can be synthesized for the 
interval plant with perturbations in that e - box. This can be done for any e < 6.321. 
The Kharitonov polynomials of the numerator are 

K n i = K n 2 = a + (1 — e) 

= K n4 = a + ( 1 + c) 


The plants corresponding to the E-segments are 

^ and 4=-; i = l,3;j = 1, 2,3,4 
where Sjj are the four segments 

[K dl (a),K d2 (a)}, [K dl (a), K^a)}, [K d2 (a),K d4 (a)], [K^a), K d4 {a)\ 

and S n is the segment [iiT„i(a), if n3 (s)]. 

By segment [Kj 1 (s), Kj 2 (s)], is meant all convex combinations of the form 

(l-A)K dl (s) + A^ 2 (s), A < 1. 


Here, it can be seen that there are 12 plants which are to be searched for the maximum 
magnitude. 

Let e = 3. MATLAB functions are used to search through the E-segment plants 
to get the plot o maxp g p |AP(ju>)| vs w as shown in Figure 2. This quantity is found 
to be 0.189. Therefore r = 0.189 Choosing u = 0.3 < 1, we have 


Therefore, 


9 = 


u 

T 


0.3 

0.189 


1.587 


C(a) 


9 _ **(*) 

1 — P 0 q r(a) - P 0 (*Ms) 
1.587(s 3 + 8s 2 + 22s + 20) 
s 3 + 8s 2 + 20.413s + 18.413' 
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Figure 2: Maximum Unstructured Perturbation vs frequency for e = 3 

So, for the stable plant case, the controller can be parametrized in terms of an arbi- 
trary Schur function u(s). This C(a) stabilizes the interval plant P(s) with coefficients 
varying in the intervals 

a G [-2,4] 
b G [5,11] 
c G [19,25] 
d G [17,23]. 

This can be verified by checking if C(s) stabilizes the plants corresponding to the 
E-segments. The numerator and denominator Kharitonov polynomials are: 

K nl = K n2 =s- 2 
K n a = K n 4 = 3 + 4 
Kji — s 3 + 11s 2 + 19s + 17 
K d2 = s 3 + 11s 2 + 25s + 17 
K d 3 = s 3 + 5s 2 + 19s + 23 
K d 4 = s 3 + 5s 2 + 25s + 23. 

The numerator and denominator Kharitonov segments are 

S n = (1 - X)K nl + XKm = s + (6A - 2) 

S d i = (1 — X)K d \ + \K d 2 = s 3 + 11s 2 + (6A + 19)s + 17 
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S d 2 = (1 — X)K d i + XKj 3 = s 3 4- (11 — 6A)s 2 + 19a 4- (17 4- 6A) 

S d 3 = (1 — X)Kdi + \K d A = s 3 4- (11 — 6A)a 2 + 25s + (17 + 6A) 

S d 4 = (1 - X)K d3 + \KdA = s 3 + 5s 2 + (19 + 6A)s + 23 

where A G [0, 1]. So, the E-segments are 

Pi = {K nl ,Sdi) 

P 2 = ( K n i,Sd2 ) 

P 3 = {K nl> S d a) 

Pa = (K nU S d4 ) 

Ps = (K^Sdi) 

Pa = {KjAiSdi) 

Pi = (K^Sdt) 

Ps = (K n3 ,Sd4) 

P 9 — (S n ,Sdi) 

P 10 = {S n) Sd2) 

P 11 = (S n) Sds) 

P 12 = (S n , Sd 4 ) ■ 

A MATLAB function was used to check if the controller C^s) stabilizes these seg- 
ments. C(s ) was found to stabilize all the segments and hence, will robustly stabilize 
the system for e = 3. 

For e = 6, the plot of ma xp e7 > |AP(ju;)| vs w is shown in Figure 3. 


max I AP(?'w)| 
Per v 


1.453 


For u = 0.3, 


Then 


9 = 


0.3 


= 0.207. 


C(s) = 


1.453 

0.207(a 3 -I- 8s 2 + 22s + 20) 


1 - P 0 q a 3 + 8s 2 + 21.793s 4- 19.793 
This compensator C(s) stabilizes P(s) with coefficient intervals 


a G [-5,7] 
b G [2,14] 
c G [16,28] 
d G [14,26] 


This is again verified by checking if (7(a) stabilizes another set of E- segments which 
correspond to the above coefficient intervals. 


So, it can be seen from the above example that by choosing a Schur u(s) it is 
always possible to come up with a robustly stabilizing compensator for an interval 
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Figure 3: Maximum Unstructured Perturbation vs frequency for e = 6 


plant built around a stable nominal plant as long as the coefficient perturbation, e is 
less than some f msxi so that there is no unstable pole over the whole family. 

Now let us consider an example for an unstable nominal plant. 

Example 2. 


P.M 


The interval plant is 


with the intervals as 


5s + 4 5s + 4 

(s — 3)(s + 5) s 2 + 2s — 15 


pm 


5s + a 
s 2 + bs + c 


a G [4 — e, 4 + e] 
b G [2 — e, 2 + e] 
c G [— 15 — e, — 15 + e]. 

The Kharitonov polynomials for the numerator and denominator are: 


K n i 

= K n2 

= 5s 

+ ( 4 - 6 ) 

K n , 

= K n 4 

= 5s 

+ (4 + 6) 

K d i 

= S 2 + 

( 2 - 

e)s 

— (15 + e) 

K d2 

= S 2 + 

(2 + 

e)s 

- (15 + e) 

K d 3 

= a 2 + 

(2- 

e)s 

+ (-15 + £) 

K d 4 

= s 2 + 

(2 + 

e)s 

+ (-15 + «) 
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Using Routh-Hurwitz criterion, it can be shown that the number of unstable zeros 
of Kji, i = 1) 2, 3, 4 (and hence the number of unstable poles of the family of plants) 
does not change for any e < 15. 

Now,the Blaschke product, 


B(s) 


3-s 
3 + s 


h(>) 




5s + 4 

(s + 3)(s + 5) 


5.3 + 4 
(3 + 3)(3 + 5) 


-0.395 


We have to find a Schur function u(s) such that 


u(3) = 


K3) 

^o(3) 


— r (3) 

0.395 


and 

|u(3)| < 1 => |r(3)| < 0.395 

Here, if the design is to be done with a constant r, it is necessary that r < 0.395 
for a robust stabilizer to exist. Therefore r max = 0.395. A plot of r vs e is shown in 
Figure 4. From the plot, € max comes out to be equal to 2.83. So the maximum e for 
which the plant can be robustly stabilized, with a constant r is 2.83. Let e = 2.83. 
Then, from the plot of ma xp G ? \AP(ju)\ vs lj } as shown in Figure 5, we have 


max |AP(?u;)| 
Pzv ' 71 


= 0.395 = r. 


Hence u(3) = 1 = u(s). Now 

u 

q= - 

T 


1 

0.395 


-2.532 


q(s) = B(s)q(s) = 


2.532 (s - 3) 
(s + 3) 


Therefore, 

r(a s g( J ) 

L[ * } ~ 1 - P 0 (s)q(s) 

This £7(s) stabilizes P(s) with 


2.532(s + 5) 
(s - 1.584) ‘ 


a G [1.17,6.83] 
b E [-0.83,4.83] 
c E [-17.83,-12.17], 
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This is tested by checking if (7(s) stabilizes the E-segments of P(s). The Kharitonov 
polynomials are: 

K n , = K n2 = 5a + 1.17 
K „3 = K n 4 — 5s 4- 6.83 
K dl = s 2 - 0.83s - 17.83 
K d 2 = s 2 + 4.83s - 17.83 
Kjz = s 2 - 0.83s - 12.17 
K d 4 = s 2 + 4.83s- 12.17 

and the Kharitonov segments are: 

S n = 5s + (5.66 A + 1.17) 

S dl = a 2 + (5.66A - 0.83)s - 17.83 
S d 2 = s 2 - 0.83s -f (5.66A - 17.83) 

S d 3 = s 2 + 4.83s + (5.66A - 17.83) 

S d 4 = s 2 + (5.66A - 0.83)s - 12.17. 

So the E-segments are: 

Pi = (K nl ,S di ) 

Pa = (K nl ,S d2 ) 

Pa = {Kni 1 S,a) 

Pa — {K n \ > S<ia) 

Pb — (Kn3,S d i) 

Pe = (Kn3, S d2 ) 

P 7 = (K^Sa) 

Pa = (K^,S d 4 ) 

Pa = (S n , K d i) 

P\ 0 = (SmKito) 

P 11 = {S n , K<a) 

P\ 2 = ( S n ,K d 4 ) 


Again MATLAB is used to verify that C7(s) stabilizes all the segments. 
For e = 4: 

11 11 

= 0.578. 


Since 


max|AP(ju/)| 


^l AP 0’ w )l 


> 0.395, 


hence r — 0.578 is not the right choice. With the plot of ma xp £ p |AP(ja;)| as in 
Figure 6 in mind, r(s) is loop - shaped to approximate it from above. One such r(s) 

lb ' 2.52 ( s + 0.6) 

^ ~ (s + 1.3)(s + 2.4)' 
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Table 2: Fenyves Array 


3 oo 


-0.987 0 

u(s) 

0.987 

u 2 (s) 


Figure 7 shows the plot of r(s). Then 


r(3) -2.52 x 3.6 

' ” p 0 (3) “ 4.3 x 5.4 x 0.395 


= -0.987. 


Since r(s) is of relative degree 1, another interpolation condition on u(s) is u(oo) = 0 
so that a proper controller is obtained. 

Now, the second row of the Fenyves array will be 


■u(oo) + 0.987 s + 3 
2 1 + 0.987u(oo) .s — 3 


= 0.987. 


Hence the Fenyves array is as shown in Table 2. 

Since the elements of the Fenyves array are of modulo less than 1, hence an SBR 
function u(s) exists which interpolates to 

u(3) = —0.987, u(oo) = 0 

This can also be verified by checking for the non-negative definiteness of the Pick 
matrix. Here 


a x = 3, 
/3a = -0.987, 


Q£j = OO 
02=0 


so that the Pick matrix is 


P = 


1-0.987 2 

6 

0 


0 

0 


4.305 x 10" 3 0 
0 0 


Since P is symmetric and has non-negative eigenvalues, hence it is non-negative def- 
inite. This also confirms the existence of u(s). 

Now, u 2 (s) can be parametrized in terms of an arbitrary Schur function u 3 (s), 

i.e., 

0.987 — u 3 (s) 


u 2 (s) = 


0.987 + u 3 (s) ( 

1— OO j 
k -+oo ) 

1 + 0.987 | 

(m or/ 
B — OO 

k*+°° - 

) u a(s) 


Choosing u 3 (s) = 0, we get = 0.987. 


15 


Then 


u(a) 




9 ( s ) 


-0.987 + 0.987 fef) -229.535 

1 + 0.987 ( s — ooa + 00 ) 1 * 3 ( 5 ) a + 229.574 

u(a) —91.0853(5 + 1.3)(« + 2.4) 
r(s) " (a + 0.6)(a + 229.574) 

_ 91.0853(5 + 1.3)(a + 2.4)(a — 3) 

} {a + 0.6)(« + 3 )(a + 229.574) ' 


Therefore, the robust stabilizer is 


C(a) 


9(f) 

l - P 0 («)g(a) 


91.0853(5 + 1.3)(j + 2.4)(a + 5) 
(a - 217.47 )(a + 2.687)(5 + 0.53) 


for the interval plant 


P(a) 


5a + a 
a 2 + ba + c 


a G [0, 8] 
b G [-2,6] 
c G [-19,-11]. 


The stability check is again verified with the E-segments. In this way, e can be 
increased in steps and for each e an attempt is made to find an u(s) which is Schur 
and satisfies the interpolation constraints. Then a robust stabilizer can be obtained 
from this u(s). 


In the case when the nominal plant is non-minimum phase, the magnitude of 
P 0 (oi) is lower, i.e. |u(ai)| is higher for a Po(^) with zeroes in the LHP than for 
a Pq(s) with zeroes of the same magnitude but in the RHP. So r(a») needs to be 
lower in the non-minimum phase case to keep |u(o£i)| less than 1, i.e. less amount of 
unstructured perturbation can be taken care of. 

Let us consider Example 2 again with the zero of the nominal plant in the RHP. 

Example S. 

5s -4 

s 2 + 2s - 15 
5 s + a 

s 2 + bs + c 

where 


Po («) = 
p(.) = 


a £ [—4 — e, — 4 + e] 
b £ [2 — e, 2 + e] 

c £ [ — 15 — c, — 15 + c] 
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The Kharitonov polynomials for the denominator remain the same but those for the 
numerator are: 


K n i = K n2 = 5s + (—4 — e) 
Kr* = K n4 = 5s + (—4 + e) 


The initial bound on e so that the number of RHP poles of the plant does not change, 
is also the same, i.e. 15. 


B(s) 

Po(s) 


3 -s 

3 + s 
—5 s + 4 

(s + 3)(s + 5) 


so that P 0 (3) = —0.23. So we find that |P 0 (3)| has decreased due to the zero of the 
plant being in the RHP. Therefore 


u(3) 


-r(3) 

0.23 


and it is required that 


r(3)| < 0.23 


for a robust stabilizer to exist, i.e., r mKX = 0.23. Again a plot of r vs. e, as shown in 
Figure 8, can be made to obtain e max . It comes out to be 1.8. Here, for e = 2, 


max|AP(ju;)| 


= 0.251 > 0.23 


Hence r = 0.251 is not the right choice. So poles and zeros are chosen so that r(s) 
approximates maxp e p |AP(ju;)| from above, 


1.591(s + 0.872) 
^ “ (s + 1.387)(s + 5)" 


This is shown in Figure 9. Then 


u(3) = 


~ r(3) 

0.23 


-0.763. 


The other interpolation condition is u(oo) = 0. So the Fenyves array is as shown in 
Table 3. The elements of the Fenyves array are of modulo less than 1 which guarantees 
the existence of an u(j). 

Here the Pick matrix comes out to be 


P = 


0.0696 

0 


0 

0 


which is obviously non-negative definite as it has non-negative eigenvalues. So u(s) 
exists. Choosing the free parameter function u 3 (a) as 0, we have u 2 (s) = 0.763. Then 


-0.763 + 0.763 fef ) -10.952 

1 - 0.763 (f^f) 0.763 _ * + H-354 
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Table 3: Fenyves Array 


3 oo 


-0.763 0 

u(s) 

0.763 

u 2 (s) 


From this the robust stabilizer is found out to be 

6.884(s + 1.387)(s + 5) 

^ ~ (s + 1.287 )(s - 17.461) 

which is found to stabilize all the E-segments of the given interval plant with coefficient 
intervals 


a G [-6,-2] 
b G [0,4] 
c G [-17,-13] 


The next example is given to show the element of conservativeness involved in 
this method. The choice of constant r(s) can limit e max to pretty small values in some 
cases. 


Example 4 • 

PoM 

The interval plant is 
with the intervals as 


30s -f 10 30s -|- 10 

(s + 2)(s - 3)(s - 2) = s 3 — 3s J — 4s + 12 

PM = 30 »+° 

s 3 + bs 7 + C3 + d 


cl G [10 — c, 10 -|- e] 
b G [—3 — e, —3 + e] 
c G [-4 - e, -4 + e] 
d G [12 - e, 12 + c]. 

The Kharitonov polynomials for the numerator and denominator are: 

K nl = K n2 = 30s + (10 - e) 

K n3 = K n4 = 30s + (10 + e) 

K d i = s 2 + (2 - e)s - (15 + e) 

K d2 = s 2 + (2 + e)s-(15 + e) 

Kjs = s 2 + (2 - e)s + (-15 + e) 

K d 4 = s 2 + (2 + e)s + (-15 + e). 



19 



Table 4: Fenyves Array 


3 

OO 


r 

0.875 

r 

0.667 

u(s) 


1.0395r 

u 2 (a) 


0.584— r 2 


Using Routh-Hurwitz criterion, it can be shown that the number of unstable zeros 
of Kdi, i = 1, 2, 3, 4 (and hence the number of unstable poles of the family of plants) 
does not change for any e < 6. 

Now, the Blaschke product, 


*(•)- 


Po(s) = B(s)P 0 (s) = 


(s — 3 )(j — 2) 

(s + 2 ){a + 3) 

30s + 10 


(s + 2) 2 (s + 3) 


so that Po(2) = 0.875 and P 0 (3) = 0.667. We have to find a Schur function u(s) such 
that 


u(2) = 

u(3) = 

for a constant r. It is necessary that 


r 

PJ2) 

r 

M 3 ) 


T 

0.875 

T 

0.667 


|u| < 1 =» |r| < 0.667. 


Now the Nevanlinna-Pick condition needs to be applied to get r max . First, the Fenyves 
array is formed, each element of which is a function of r. The element of the second 
row is given by 


u 2 (3) 


0.667 


0.875 


1 — 

A 0.584 

1.0395r 


r3 + 2i 
.3-2. 


0.584 -r 2 ' 


So the Fenyves array is 

Now T’ma* is obtained by equating the modulus of the element in the second row 
of the Fenyves array to 1, and solving the resulting quadratic equation. 

1.0395r ^ 

0.584 - r 2 ~ 


The above equation when solved under the constraint of |r| being already less than 
0.667, gives t < 0.404. So we have 7* max = 0.404 and u 2 (3) = 1. Now using Nevanlinna 
Interpolation mapping functions, we get 


u(s) 


a - 0.736 
a + 0.736 
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Figure 10: r vs. e 


from which, we get 

2.475(s — 0.736)(s + 2) 

^ ~ (* + 12.366)(a + 0.365)’ 

From the e - r plot as shown in Figure 10, we find that e max = 0.85 which is pretty 
small. Coefficient perturbations greater than this can be stabilized only if the choice 
of r(s) is less conservative. Figure 11 shows the maximum unstructured perturbation 
vs frequency plot for e = 0.85. This C'(a) stabilizes P(s) with 

a € [9.15,10.85] 
b G [-3.85,-2.15] 
c 6 [-4.85,-3.15] 
d € [11.15,12.85] 

This is again tested by checking if C(a) stabilizes the E-segments of P(s).The test 
was successful, i.e., C(s) indeed stabilized all the E-segments. 

So, we find that the Nevanlinna-Pick synthesis method developed in [2] can be 
used to form a robust stabilizability condition in the parametric framework. Though 
this condition is dependent on the procedure and is not absolute, it can be used as 
an important synthesis tool for interval plants within reasonable limits. 
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Figure 11: Maximum Unstructured Perturbation vs frequency for e — 0.85 

5. CONCLUSION 

The methods used in this paper for establishing robust stabilizability conditions for 
interval plants and the synthesis of robust stabilizers are well-known in the area of 
H-Infinity Robust Control. The synthesis problem in the Parametric Robust Control 
area has never been attacked before. So the parametric problems are modified and 
adapted to a form so that the H-infinity methods could be easily used for them. The 
results obtained are conservative and procedure - dependent. Inspite of all this, the 
methods developed in this paper are neat and simple and shows one approach to 
solving the synthesis problem. 
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